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ABSTRACT 


SNPs from nrDNA and ep trnC-trnD were analyzed for J. 
excelsa, J. polycarpos var. polycarpos, J. p. var. seravschanica and J. 
p. var. turcomanica and compared to terpene and RAPDs data. These 
data, taken together, support the continued recognition of J. excelsa 
and J. polycarpos as separate species as well as the three varieties of J. 
polycarpos: var. polycarpos, var. seravschanica and var. turcomanica. 
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The taxonomy of J. polycarpos K. Koch. from Armenia, 
Kazakhstan, Pakistan and Turkmenistan has been examined using leaf 
oil compositions and DNA fingerprinting (Adams, 1999; Adams, 
2001). The compositions of the volatile leaf oils (Adams, 2001) are 
given in Table 1. Notice the large amounts of a-pinene in the plants 
from Armenia, Kazakhstan and Turkmenistan. Myrcene is a large 
component in the oils from Kazakhstan and Pakistan (J. p. var. 
seravschanica). Several compounds distinguish J. excelsa: decadienal 


Phytologia (August 2008) 90(2) 209 


isomer (KI 1312), trans-cadina-1(6),4-diene, cubebol, 1-epi-cubenol, 
and KI 1666 (Table 1). Compounds that distinguish J. polvcarpos 
(including J. p. var. turcomanica and J. p. var. seravschanica for this 
discussion) are: hexyl 3-methyl butanoate, à-elemene, y-cadinene, 
elemol, germacrene B, germacrene D-4-ol, a & f-eudesmols and KI 
1688 (Table 1). Several diterpenes are unique to J. procera (Table 1) 
and show its separation from J. excelsa and the other junipers. 


The trend in the volatile leaf oils is seen in figure 1. The leaf 
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Figure |. PCO of J. excelsa, J. polycarpos var. polycarpos, J. p. vat. 
seravschanica, J. p. var. turcomanica and J. procera based on 106 leaf 
terpenoids scored as present (+) or absent (-). Based on Adams (2001) 
data. 
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terpenoids clearly separate J. excelsa and J. procera from J. 
polycarpos. The varieties of J. polycarpos are not well resolved. This 
is seen in the raw data in table 1. 

PCO analysis based on RAPDs (Adams, 2001) shows (Fig. 2) 
a very similar pattern to that seen with terpenoids. Notice that again, J. 
excelsa and J. procera are very well resolved from J. polycarpos (Fig. 
2). There appears to be more slightly more separation of the J. 
polycarpos varieties in the RAPDs data (Fig. 2) than in the terpenoids 
(Fig. 1) but the overall trend is very similar. 
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Figure 2. PCO based on 126 RAPDs. Notice that 43% of the variance 
separated J. excelsa/J. procera from J. polycarpos/turcomanical 
seravschanica (axis 1) and 15% of the variance separates J. procera 
from all other taxa on axis 2. Based on Adams (2001) data. 
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The pattern is simplified by removing J. procera from the 
RAPDS data, as seen in figure 3. Juniperus excelsa is still well 
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Figure 3. PCO of J. excelsa, J. polycarpos, J. p. var. seravschanica and 
J. turcomanica using 106 RAPDs. Based on Adams (2001) data. 


resolved. However, there is now some separation between J. 
polycarpos, J. p. var. seravschanica and J. tircomanica (Fig. 3). 


Removing J. excelsa from the data set and re-analyzing the 
RAPDs data gave a clearer picture of the pattern among the J. 
polvcarpos varieties (or populations). This resulted in four groups, 
each well resolved. The J. p. var. seravschanica populations are 
resolved into the Pakistan and Kazakhstan sites. 


Should these entities be recognized as varieties of J. 
polvcarpos or do they merely represent geographical interspecific 
variation? 
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Figure 4. PCO of four taxa of J. polycarpos based on RAPDs. Based 
on Adams’ (2001) data. 


Farjon (2005, p. 291) treated J. polvcarpos as a variety of J. 
excelsa (J. excelsa var. polycarpos (K. Koch) Takht.) and treated J. p. 
var. seravschanica and J. p. var. turcomanica as synonyms of J. 
excelsa var. polvcarpos. Clearly, neither the terpenoids nor RAPDs 
support Farjon's merging of J. polvcarpos and J. excelsa. 


Farjon (2005, p. 343) states his philosophy as "l consider 
species based on the chemistry of terpenes and/or RAPD analysis as 
based on inconclusive evidence" although he allows that "DNA 
sequence data certainly can (which is the main reason for their 
'superiority)" (Farjon, 2005, p. 232). 


The goal of sequencing a single gene is now easily 
accomplished but these data are proving to be more difficult to 
interpret than perhaps imagined. For example, Syring et al. (2007) 
examined Pinus species utilizing sequences from three nuclear genes 
(AGP6, cesAl1, LEA-like). They found that none of the three genes, 
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analyzed separately, placed the multiple accessions of P. strobus, P. 
monticola and P. chiapensis into monophyletic clades. We suggest 
that additional data from multiple genome sites such as RAPDs as well 
as phenotypic data such as leaf terpenoids can help to complement 
single gene sequences. In cases where sequence data alone results in 
multiple alternative scenarios, those additional data might be able to 
help discriminate between them. 


Recent DNA sequence phylogenetic research of Juniperus 
(Schwarzbach et al, in prep.) has shown (Fig. 5) that J. polycarpos 
(Armenia, cf. J. p. var. polycarpos, above) is 100% supported as being 
separate from J. procera and J. excelsa and not forming a 
monophyletic group. The thesis of Farjon (2005) that J. polycarpos is 
a variety of J. excelsa is not supported by these data. In fact, the 
sequencing data is in full agreement with both the terpenoid and RAPD 
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Figure 5. Bayesian tree based on nrDNA and cp trnC-trnD sequences. 
Juniperus polycarpos is resolved from J. procera, J. excelsa, J. 
foetidissima and J. thurifera. Numbers above the branch points are 
posterior probabilities on a percent basis. 
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data (above) that show J. excelsa, J. polycarpos, and J. procera being 
well differentiated. In fairness to Farjon (2005), it should be conceded 
that if only morphological data are utilized, then one could readily 
make a case for the merging of J. excelsa and J. polvcarpos (and 
perhaps J. procera). However, morphology can be misleading when 
used solely for species circumscriptions, as has been shown in several 
other studies. Cryptic speciation has been found in J. deltoides R. P. 
Adams (Adams et al. 2005) and J. maritima R. P. Adams (Adams 
2007), to name but two cryptic juniper species. 


Schwarzbach et al. (in prep) have analyzed one accession per 
species for J. excelsa, J. polycarpos, and J. procera. As a result, the 
relationships of the species were established in a basic framework. 
However, the sampling in this previous study did not allow any 
assessments of intraspecific variation or the monophyly of the taxa 
involved. The purpose of the present paper is to re-examine the 
taxonomy of J. excelsa, J. polycarpos and its varieties using SNPs 
from sequence data (nrDNA and cpDNA tmC-trnD) by adding 
multiple accessions for each taxon. 


MATERIALS AND METHODS 


Specimens used in this study: J. exce/sa, Adams , 8785-8786 - 
7 km w of Lemos, Greece; Adams 9433-9435, 40 km n of Eskisehir, 
Turkey; J. p. var. seravschanica, Adams 8224-8226, 2 km s 
Dzhabagly, Kazakhstan (not Kyrgystan as previously reported, Adams, 
1999); J. p. var. seravschanica, Adams 8483-8486, Quetta, 
Balochistan, Pakistan; J. turcomanica, Adams 8757-8760, Kopet Mts., 
Turkmenistan. Voucher specimens for all collections are deposited at 
BAYLU. 

One gram (fresh weight) of the foliage was placed in 20 g of 
activated silica gel and transported to the lab, thence stored at -20° C 
until the DNA was extracted. DNA was extracted from juniper leaves 
by use of a Qiagen mini-plant kit as per manufacturer's instructions. 


SNPs obtained from DNA sequencing 

ITS (nrDNA) and trnC-tmD amplifications were performed in 
30 ul reactions using 6 ng of genomic DNA, 1.5 units Epi-Centre Fail- 
Safe Taq polymerase, 15 ul 2x buffer E or K (final concentration: 50 
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mM KCl, 50 mM Tris-HCl (pH 8.3), 200 uM each dNTP, plus Epi- 
Centre proprietary enhancers with 1.5 - 3.5 mM MgCl, according to the 
buffer used) 1.8 pM each primer. All 12 (A-L) of the Epi-Centre's 
buffers were screened and buffer K gave the cleanest, most-abundant 
amplification for both ITSA/ITSB and buffer E was best for trnC-trnD 
(CDIOF/CD3R) primers. However, buffers D, F, G, H, and J were 
nearly as good as buffer E or K. 

Primers (5-3): for nrDNA: ITSA and ITSB primers from 
Blattner (1999), for trnC-tmD: CDIOF and CD3R, see Adams et al. 
(2008) The following PCR conditions were used: MJ Research 
Programmable Thermal Cycler, 30 cycles, 94°C (1 min.), 50°C (2 min.), 
72°C (2 min.), with a final step of 72°C (5 min.). The PCR reaction 
was subjected to purification by agarose gel electrophoresis (1.5% 
agarose, 70 v, 55 min.). The nrDNA primers (ITSA, ITSB) produced a 
band of approx. 1120 bp. The internal trnC-trnD primers, CD10F- 
CD3R produced a band of approx. 850 bp. In each case, the band was 
excised and purified using a Qiagen QlAquick gel extraction kit. The 
gel purified DNA band with the appropriate primer was sent to McLab 
Inc. for sequencing. Sequences for both strands were edited and a 
consensus sequence was produced using Chromas, version 2.31 
(Technelysium Pty Ltd.). Alignments were made using Clustal W and 
then manually corrected. Indels were coded with a "-" for the first 
nucleotide and "I" for succeeding nucleotides such that an indel was 
treated as a single mutation event. Overall sequences have been 
deposited in GenBank (Schwarzbach et al., in prep.). 


SNPs analyses 

Aligned data sets (nrDNA and trnC-trnD) were analyzed by 
CLEANDNA (Fortran, R. P. Adams) to remove invariant data. 
Mutational differences were computed by comparing all SNPs, divided 
by the number of comparisons over all taxa (= Gower metric, Gower, 
1971; Adams, 1975). Principal coordinate analysis was performed by 
factoring the associational matrix using the formulation of Gower 
(1966) and Veldman (1967). A minimum spanning network was 
constructed by selecting the nearest neighbor for each taxon from the 
pair-wise similarity matrix, then connecting those nearest neighbors as 
nodes in the network (Adams, 2004). 
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RESULTS AND DISCUSSION 


Sequencing nrDNA (ITS region) resulted in 1210 bp of data. 
Aligning sequences for J. excelsa (5 individuals, Greece, Turkey), J. 
polycarpos var. polycarpos (4 individuals, Armenia); J. p. var. 
seravschanica (4 individuals, Kazakhstan, 4 individuals, Pakistan) and 
J. p. var. turcomanica (4 individuals, Turkmenistan) revealed 14 SNPs 
among these individuals. PCO of these individuals gave three 
significant eigenroots accounting for 74.4%, 9.1% and 6.7% of the 
variance among individuals. PCO ordination shows (Fig. 6) two major 
groups to be present: J. excelsa - J. p. var. turcomanica and J. p. var. 
polvcarpos and var. seravschanica. 
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Figure 6. PCO based on 14 nrDNA SNPs. Dashed lines show the 
minimum linkage between groups. Numbers above the dashed lines 
are the number of SNP events separating the groups. Equally spaced 
lines denote identical DNA sequences. 
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Sequencing the partial trnC-trnD sequence resulted in 877 bp 
of data when utilizing the same genomic DNA as above. Aligning 
these 21 sequences revealed 14 SNPs. Factoring the association matrix 
resulted in 3 eigenroots that accounted for 49.1%, 43.9% and 3.3% of 
the variance among these individuals in their partial trnC-trnD SNPs. 
Ordination of these individuals reveals (Fig. 7) three groups: J. excelsa; 
J. polycarpos - J. p. var. turcomanica; and J. seravschanica. These 
groups are separated by 7 and 8 SNPs. Notice (Fig. 7) that there was 
no variation among individuals of any taxon except J. excelsa. 
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Figure 7. PCO ordination based on 14 SNPs from tmC-tmD sequences. 
Dashed lines show the minimum linkage between groups. Numbers 
above the dashed lines are the number of SNP events separating the 
groups. Equally spaced lines denote identical DNA sequences. 


Combining the 14 nrDNA SNPs and 14 tmC-tmD SNPs for a 
PCO analysis resulted in 3 eigenroots of 52%, 30% and 12%. 
Ordination (Fig. 8) shows 4 well-defined groups, each separated by 7 
to 9 SNPs. It is interesting that the 8 samples of J. p. var. 
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seravschanica (4 from Kazakhstan and 4 from Pakistan) had identical 
sequences for both nrDNA (1210 bp) and cp trnC-tmD (877 bp). 
Juniperus p. var. polycarpos and J. p. var. turcomanica both had a 
single SNP within their samples. In contrast, J. excelsa had 2 SNPs 
within its samples. 
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Figure 8. PCO using SNPs from both nrDNA and trnC-trnD. Note the 
four well defined groups. Dashed lines show the minimum linkage 
between groups. Numbers above the dashed lines are the number of 
SNP events separating the groups.  Equally spaced lines denote 
identical DNA sequences. 


Graphic summaries of morphology, terpenes (+/- basis), 
RAPDs, nrDNA, trnC-trnD and combined nrDNA + trnC-trnD reveal 
(Fig. 9) general agreement between morphology, terpenes, RAPDs and 
combined nrDNA + tmC-tmD classifications. The combined nrDNA + 
trnC-trnD SNPs showed the largest differences between J. polycarpos 
var. polycarpos, J. p. var. seravschanica and J. p. var. turcomanica of 
any data set. Using only nrDNA or tmC-trnD SNPs would lead to very 
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different taxonomies in this study. The concordance of terpenes, 
RAPDs and combined nrDNA + trnC-tmD classifications seems to 
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Figure 9. Graphic summaries of morphology, terpenes, RAPDs, 
nrDNA, trnC-trnD, and combined nrDNA + tmC-trnD data. E = J. 
excelsa, P = J. p. var. polycarpos, T = J. p. var. turcomanica, Sk = J. p. 
var. seravschanica, Kazakhstan, Sp = J. p. var. seravschanica, 
Pakistan. 
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provide the strongest evidence that J. polycarpos is composed of three 
genetically distinct (but scarcely distinct in morphology) taxa, 
supporting the continued recognition of these taxa at the variety level. 
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Table |. Comparisons of the percent total oil for leaf essential oils for J. 
excelsa - Greece (EG), J. excelsa, Tbilisi Botanic Garden (ET), J. polycarpos 
var. polycarpos: Armenia, L. Sevan (AS) J. p. var. turcomanica, 
Turkmenistan, Kopet Mts. (TK), Alma Ata Botanic Garden (ex. Ashgabad, 
Turkmenistan, TA); J. p. var. seravschanica: Kazakhstan, Talasskiy Mtns. 
(KT), Pakistan, Quetta (PQ), and J. procera, east Africa (PR). Components that 
tend to separate the species are highlighted in boldface. From Adams (2001). 


KI Compound EG ET AS TK TA KT PQ PR 
926 tricyclene 0.1 0.1 01 O1 02 0.3 O. t 
931 u-thujene - - t t t 0.6 0.4 t 
939 a-pinene 22.5 26.5 68.4 68.8 59.7 44.4 15.5 12.5 
953 a-fenchene 0.2 t t t t - 020l 
953 camphene 05 10 02 O01 05 0.5 02 0.1 
957 thuja-2,4(10)-diene 0.1 - t t t - - - 
975 verbenene t - - - - - - - 
976 sabinene t 01 02 01 04 0.9 0.5 t 
978 1-octen-3-ol - - - - - - - 03 
980 B-pinene 0:6 10) 05 06 TS 2 2 TINO MEMO: 
99] myrcene 19 22 12 15 37 192 OIE 
1005 a-phellandrene 0.1 0.1 - - t 0.1 0.1 - 
1011 ô-3-carene 2.3 04 t t t > SS (ull 
1018 a-terpinene 0.1 0.1 t t t 0.1 0.1 t 
1026 p-cymene 04 02 0.1 O01 02 01 07 t 
1028 sylvestrene - - - - - - - 01 
1031 limonene 22.6 55 12 15 18 44 90 02 
1031 f-phellandrene t - - - Ol 0.5 10 08 
1032 1,8-cineole - t - - - - - t 
1050 (E)-B-ocimene t - t t t t 0.2 t 
1062 y-terpinene 0.6 0.9 02 03 06 1.4. 13 t 
1068 cis-sabinene hydrate - - - - - 0.1 0.2 - 
1068 fenchone - - t t t t t - 
1088 terpinolene 09 || 04 O05 13 rey lla Al 
1097 trans-sabinene hydrate - - - t - - - 
1098 linalool - 01 OIl t t 0.5 07 0.5 
1103 isopentyl-isovalerate - - t t t - - - 
1110 1,3,8-p-menthatriene - - - - - - - t 
1112 endo-fenchol 0.2 - - - t - - - 
1114 trans-thujone - - t t 02 0.1 - - 
1121 cis-p-menth-2-en-1-0] 0.1 - - - - t - t 
1125 chrysanthenone - - t t t - - - 


1125 a-campholenal Gi $92 02 i @2 t Ol t 
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bornyl acetate 
linalool oxide acetate 
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trans-sabinyl acetate 
decadienal isomer? 
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(-copaene 
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Kl Compound EG ET AS TK TA KT PQ PR 


1418 B-cedrene 0.9 0.5 - - - 01 02 - 
1429 cis-thujopsene 03 02 02 - - 02 04 - 
1446 cis-muurola-3,5-diene 0.2 0.6 - 0.2 - t - - 
1454 a-humulene (Q2 (2 s - - EN NI 
1458 E--farnesene 02 Ol Ol - - - 0.1 - 
1461 cis-muurola-4(14),5- 

diene - - 01 02 - - 01 - 
1466 D-acoradiene 0.1 t 0.1 - - - t - 
1473 trans-cadina-1(6),4- 

diene 0.4 0.8 - - - - - - 
1477 y-muurolene - t - 02 0.1 t - - 
1480 germacrene D 09 17 02 08 08 OI 02 03 
1491 trans-murrola-4(14),5- 

diene 04 14 - 01 QI t - - 
1493 epi-cubebol - ee OE t - 
1499 y-amorphene - - - - - - t - 
1499 o-muurolene 02 0.1 - 03 02 02 02 - 
1499 bicyclogermacrene - - t - - - t 
1502 cuparene - - - - - - t - 
1503 germacrene A - - - O1 01 - - 
1509 -bisabolene - - t - - - t 
1512 a-alaskene 0.3 t 0.1 - - - 04 - 
1513 y-cadinene - - 02 11 09 04 04 - 
1513 cubebol 0.8 26 - - - - - - 
1521 cis-calamenene t - - - - - - - 
1524 ó-cadinene Oy des uS Jub 12 i) (Quy = 
1532 E-y-bisabolene () = OIl - - - - - 
1532 trans-cadina-1(2),4- 

diene t 02. - t - t - - 
1538 a-cadinene - - - 02 02 01 0.1 - 
1549 elemol - - t 02 19 05 05 43 
1556 germacrene B - - 07 16 18 04 0.6 - 
1574 germacrene D-4-ol - t 04 30 15 05 2.0 0.1 
1581 caryophyllene oxide - - - 
1587 sesquiterpene, FW220? 1.9 1.7 0.8 - - 1.0 1.7 - 
1596 cedrol 28.1 30.8 19.0 t - 146 264 - 
1606 humulene epoxide ll t - - - - - - 05 
1606 B-oplopenone t - t 
1607 4E-tridec-6-yne* - 06 - 
1627 l-epi-cubenol 1.6 2.2 t - - - - - 
1630 y-eudesmol - - - 
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KI 


1640 
1640 
1645 
1649 
1652 
1653 
1666 
1666 


1685 
1688 


1688 
1789 
1809 


1930 
1961 


1989 
2054 
2080 
2103 


2147 
2181 
2278 
2288 
2293 
2302 


2302 
2325 


Compound EG 
epi-a-cadinol t 
epi-a-muurolol t 
a-muurolol t 


B-eudesmol - 
o-eudesmol - 
a-cadinol t 
bunesol - 
unknown(57,41,85, 
79,136) 0.6 
eudesma-4(15),7-dien- 
1-f-ol - 
sesquiterpene alcohol, 
FW 222 - 
cadinol isomer - 
8-a-acetoxyelemol - 
unknown(43,79,7 1,99, 
136,252) - 
rosa-5,15-diene - 
sandaracopimara- 
8(14),15-diene t 
manoyl oxide t 
abietatriene t 
abietadiene 0.3 
diterpene, 41,79,191, 
257,F W286? - 
abieta-8(14),13(15)- 
diene* - 
diterpene,41,91, 271, 
257,F W286 - 
sempervirol - 
4-epi-abietal 0.1 
diterpene,41,55,255, 
269,FW284? - 
abieta-7,13-dien-3-one  - 
trans-totarol - 
trans-ferruginol - 


ET AS TK TA KT PQ 
03 01 07 09 OS O04 
0.5 t VEO Oil tl 
0.2 - 02 Qd t t 

- t t 0.8 t 0.1 
- t t 0.5 0.1 t 
02 02 16 32 09 O08 
- t t DA = 01 
2.4 - - - = = 
- 03 09 36 0.4 - 
- - t 1.2 - - 
0.6 - - t - - 
02 01 0.1 - o M3 
0.4 t t 0.3 - t 
22 03 07 - a | 7 
12 Oy il ET - 1.0 
0.1 - - 02 2 5 
0.3 - = = = 2 


KI = Kovat's Index on DB-5(=SE54) column. *Tentatively identified. 


Compositional values less than 0.1% are denoted as traces (t). 


Unidentified components less than 0.5% are not reported. 


